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ABSTRACT 



Context. Jets of plasma are frequently observed in the solar corona. A self-similar recurrent behavior is observed in a fraction of 
them. 

Aims. Jets are thought to be a consequence of magnetic reconnection, however, the physics involved is not fully understood. Therefore, 
we study some jet observations with unprecedented temporal and spatial resolutions. 

Methods. The extreme-ultraviolet (EUV) jets were observed by the Atmospheric Imaging Assembly (AIA) on board the Solar 
Dynamics Observatory (SDO). The Helioseismic and Magnetic Imager (HMI) on board SDO measured the vector magnetic field, 
from which we derive the magnetic flux evolution, the photospheric velocity field, and the vertical electric current evolution. The 
magnetic configuration before the jets is derived by the nonlinear force-free field (NLFFF) extrapolation. 

Results. Three EUV jets recurred in about one hour on 2010 September 17 in the following magnetic polarity of active region 1 1 106. 
We derive that the jets are above a pair of parasitic magnetic bipoles which are continuously driven by photospheric diverging flows. 
The interaction drove the build up of electric currents that we indeed observed as elongated patterns at the photospheric level. For the 
first time, the high temporal cadence of HMI allows to follow the evolution of such small currents. In the jet region, we found that 
the integrated absolute current peaks repetitively in phase with the 171 A flux evolution. The current build up and its decay are both 
fast, about 10 minutes each, and the current maximum precedes the 171 A by also about 10 minutes. Then, HMI temporal cadence is 
marginally fast enough to detect such changes. 

Conclusions. The photospheric current pattern of the jets is found associated to the quasi-separatrix layers deduced from the magnetic 
extrapolation. From previous theoretical results, the observed diverging flows are expected to build continuously such currents. We 
conclude that magnetic reconnection occurs periodically, in the current layer created between the emerging bipoles and the large scale 
active region field. It induced the observed recurrent coronal jets and the decrease of the vertical electric current magnitude. 

Key words. Magnetic fields - Sun: corona - Sun: surface magnetism - Sun: UV radiation 



1. Introduction 

Solar jets are observed at many wavelengths , such as Ha, 
UV, EUV, and X-ray s (ISchmieder et all Il988l IShibata et al. . 
1992a] ICanfield et all Il996t IChae eTail Il999t It 
20121) . They are called surges in cold 
Ha) and jets at other wavelengths 



some of the aforementioned models focus on solar flares and 
coronal bright points, the associated jet phenomenon has also 
been discussed. 



Uddi n et al 
spectral lines (e.g., 
which indicates that 



Ha surges and EUV/X-ray 
same region tend to appear re 



jets eman ating from the 



both cold (e.g., ISrivastava & Murawskil 120111: iKav shap et al., 
1201 3b and hot plasma are accelerated in the jet phenomenon. 
Magnetic reconnection is widely accepted as the central phys- 
ical mechanism. There are different magnetic reconnection 
models for coronal jets according to photospheric flow pat- 
terns and coronal magnetic topologies, for instance, the emerg- 
ing flux model (iHevvaertset all Tl 9771: IShibata et all Il992bt 
iGontikakis et all l2009t lArchontis et all 12010). the converg- 
ing flux model (Priest et al., 1994), the Quasi-Separatrix Layer 
(QSL) model dMandrini et all Il996|). and the null-point and 
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Zhans et all 


20121). The period ranges from tens of minutes to 



20081 iToroket all 
(120081) 

use emerging flux, while iPariat et all 



fan-se paratrix mod el (Moreno-Insertis et al. 
2009; Pariatet all I2009L l2010l). iMoreno-Insertis etal 



and Torok et all (120091) 

(12009! 120101) " use horizontal photospheric twisting motion to 
drive the magnetic reconnection in their simulations. Although 



several hours. It is still not clear which physical me chanism ac- 
counts for the quasi-periodicity. Pariat et ail (1201 Oh proposed a 
null-point and fan-separatrix model driven by the photospheric 
twisting motion. As a consequence of the continuous driving by 
photospheric motions, this mo del produces recurre nt jets, which 
explains the quasi-periodicity. Zhang e t alJ (1201 2l) proposed an- 
other possible mechanism that the null-point reconnection re- 
sponsible for the quasi-periodicity is modulated by trapped slow- 
mode waves along th e spine field lines, an d slow-mode wave has 
been reported in, e.g JCirtain et al.l(l2007l) . Both models consider 
magnetic topologies with null-points. But how to explain those 
events whose magnetic topologies are QSLs or bald patches? 
Moreover, the photospheric motion is believed to be the driver 
of these recurrent jets. Can we find any indicator of the quasi - 
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Fig. 1. (a)-(c) SDO/AIA EUV composite images. Each EUV image consists of three channels (red, blue, and green) in the 304 A, 
171 A, and 21 1 A channels, whose characteristic temperature corresponds to 0.05 MK, 0.6 MK, and 2 MK, respectively. The larger 
solid box marks the jet region where the 171 A fluxes are integrated. The smaller solid box on the clock indicates the position and 
field of view of the EUV image on the solar disk. The dashed box indicates the region where the background flux for the 171 A band 
is computed, (d) SDO/HM1 line-of-sight magnetic field. The black/white color represents the negative/positive polarity. The solid 
box marks the region of the foot-points of the jets, where the line-of-sight magnetic field flux is integrated. The dotted box marks 
the field of view for the analysis of the vector magnetic field. 



periodicity from the photosphere, such as from the magnetic flux 
evolution, the velocity field, or the electric current evolution? 

In this paper, we present three recurrent coronal jets that 
were observed on 2010 September 17 by the Atmospheric 
Imaging Assembly (AIA: lLemen et al1l2012l) on board the Solar 
Dynamics Observatory (SDO). In order to find a clue on the 
magnetic reconnection model and the physical mechanism ac- 
counting for the quasi-periodicity, we study the EUV flux, the 
magnetic flux, the velocity field, and the vertical electric cur- 
rent evolutions. The vector magnetic fields are observed by the 
Hehoseismic and M agnetic Imager (HMI; IScherrer et al]l2012t 
ISchou et"a"Dl2012bT) on board SDO. Section [2] presents the data 
analysis and results. Discussion and conclusion are given in 
Section[3] 

2. Data Analysis and Results 

2.1. Jets Observed in EUV 

SDO/AIA provides high cadence of 12 s, high spatial resolution 
of 1.5" (the pixel sampling is 0.6" per pixel), and high signal- 
to-noise observations nearly simultaneously in seven EUV lines, 
two UV continuums, and one white-light band. Figure Q] dis- 



plays three composite EUV images showing that three coronal 
jets recurred on the border of active region (AR) 11106. All 
the EUV images have been aligned with "aia_prep.pro" in the 
Solar SoftWare (SSW). A line-of-sight magnetic field observed 
by SDOfHMl is displayed in Figure [TjJ as an example, which 
shows that some parasitic positive polarities distribute in the 
main negative polarities. SDOfHMl has a spatial resolution of 
1" (the spatial sampling is 0.5" per pixel) and a cadence of 45 s 
for the line-of-sight magnetic field. It has also been aligned with 
the EUV images with "aia_prep.pro". 

We compute the EUV flux in the 171 A band to give a quan- 
titative representation of the evolution of the recurrent jets. The 
171 A flux is computed within a rectangle containing the full 
jet region as shown by the larger solid box (130" x 80") in 
Figure [TJi-QJ- We select two time ranges to compute the 171 A 
flux evolution, i.e., a context range from 00:00 to 10:00 UT and a 
smaller analysis range centered on the analyzed jets from 03:00 
to 05:00 UT. Due to the limitation of computation resources, the 
temporal resolution for the context range is 2 minutes and it is 
12 seconds for the analysis range. The 171 A flux evolutions for 
the two time ranges are plotted in Figures [2^ and |2j), respec- 
tively. The data gap between 06:31 to 07:13 UT is due to the 
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Fig. 2. SDO/AIA 171 A averaged and normalized flux and 
SDO/HM1 line-of-sight magnetic flux, (a) Blue stars indicate the 
171 A flux in the jet region divided by the number of pixels in 
the box (larger solid box as shown in Figure QJ-QJ) normalized 
to the background flux, which is computed in the dashed box 
shown in FigureQ] Plus signs represent the averaged background 
flux per pixel. The temporal resolution is about 2 minutes. The 
two vertical lines mark the time range for 171 A flux as shown 
in the bottom panel, (b) The temporal resolution for the 171 A 
flux is about 12 s. The three vertical lines mark the peak times of 
the 171 A flux. Dash-dotted, dotted, and solid lines indicate pos- 
itive, unsigned negative, and total unsigned line-of-sight mag- 
netic flux, respectively. 



SDO eclipse by the Earth. We have normalized the 171 A flux 
to its background level, which is computed in a quiet region as 
shown in the dashed box of Figure[]J-[3J- The 17 1 A background 
flux shows a flat curve (Figure |2^). Therefore, the variations in 
the 171 A flux curve in the jet region are mainly due to the jets 
themselves, but not the background. 

Figure [2^ shows many fluctuations or peaks in the 171 A 
flux curve during the context range, which indicates that there 
are many brightening points or jets during this time range in the 
selected region. However, the peaks during the analysis range 
are more distinct than the others. Therefore, we select this time 
range for a detailed analysis. The 171 A flux curve in Figure|2j) 
has three main peaks at 03:17 UT, 03:57 UT, and 04:22 UT, 
respectively. There are some smaller peaks before the jet at 
03:57 UT. Using the 171 A movie (attached to the online-only 
Figure |7J, we find that this jet consists of successive ejections 
during a short period. The other two events are relatively simple 
with only one ejection in each event. 

2.2. Line-of-Sight Magnetic Flux Evolution 

Figure |2j5 displays the unsigned line-of-sight magnetic flux evo- 
lution. The line-of-sight magnetic flux is integrated in a small 
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Fig. 3. (a) SDO/AIA 171 A image overlaid by the line-of-sight 
magnetic field observed at the same time as that of the 171 A 
image. Gray/black lines mark the contour levels of the posi- 
tive/negative polarities. The contour levels of B z are -1200, - 
900, -600, 300, and 600 G. (b) Map of photospheric trans- 
verse velocities derived from LCT analysis with the HMI mag- 
netograms, employing a FWHM correlation window of 5" dur- 
ing the time interval 02:25-02:50 UT. Arrows represent veloc- 
ities and the background is the line-of-sight magnetogram with 
positive (white) and negative (black) polarities. Black/white ar- 
rows are located on positive/negative polarities. North/south is 
up/down, west/east is towards right/left. The magnetic polari- 
ties involved in the jets are identified with annotations (pi and 
p2 for parasitic positive, nl and n2 for parasitic negative, while 
NO and Nl for the main following negative polarities). The cir- 
cle highlights the region of interest where the diverging flows 
are found. A movie showing the evolution of the line-of-sight 
magnetic field from 02:00 to 06:00 UT is available in the online 
journal. 



rectangle containing mainly the foot-points of the jets with a size 
of 48" x28" as shown by the solid box in FigureQJl. The box size 
for the 171 A flux computation is larger than that for the mag- 
netic flux because, on one hand, jets are events occupying areas 
in the corona that are more extended than the triggering mag- 
netic field region in the photosphere below. On the other hand, 
the parasitic polarities involved with the jets are surrounded by 
main polarities of the active region, which would bias the com- 
putation of the line-of-sight magnetic flux if the box is too large. 
The positive magnetic flux increases from 02:00 to 03:00 UT 
and evolves almost steadily from 03:00 to 05:00 UT (Figure |2j), 
dash-dotted curve). The negative flux curve is dominated by the 
main negative polarity of the active region and some of the flux 
is close to the border of the selected region (Figure [TJl), then, due 
to convective motions of the photospheric plasma, magnetic flux 
might be displaced and no longer captured by the selected area. 
For this reason, the unsigned negative flux evolution (Figure [2J3, 
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dotted curve) is not necessarily related to the jetting events. It 
evolves steadily from 02:00 to 03:00 UT and decreases during 
the recurrent jets from 03:00 to 05:00 UT. Comparing the mag- 
netic flux evolution with the 171 A flux evolution, we do not find 
a clear relationship between them. The unsigned magnetic flux 
in the negative polarity, either decreases (for the jets at 03:17 
and 04:22 UT) or increases (for the jet at 03:57 UT) before a jet. 
There is also no significant evolution of the positive magnetic 
flux related to the jets. This result is coherent with a progres- 
sive storage of magnetic energy in the corona, reaching a critical 
point, and release part of the free magnetic energy without af- 
fecting the magnetic flux crossing the photosphere. 

Figure [3}i shows that the foot-points of the EUV jets are lo- 
cated where parasitic polarities appear, as follows. As shown in 
Figure fljl, AR 11106 consists of a leading positive and a fol- 
lowing negative polarity, which has two strong concentrated po- 
larities NO and Nl. The parasitic polarities noted as pi, p2, nl, 
and n2 in Figure [3j> are present in the following negative po- 
larity of AR 11106. The positive parasitic polarities pi and p2 
separate westward from the eastern negative polarity (nl and 
n2). Polarities nl and n2 have almost merged with the follow- 
ing main polarity Nl of AR 11106, while p2 is pushed against 
another part of the following main polarity (NO). The evolution 
of these polarities can be seen in the movie of the line-of-sight 
magnetic field attached to Figure [3] (available in the online jour- 
nal). 

We further study the transverse flows of photospheric 
magnetic features w it h loca l correlation tracking techniques 
(INovember & SimonL 1 19881 hereafter LCT). By using a 
Gaussian tracking window of full width at half maximum 
(FWHM) of 5", we compute the proper motions of magnetic ele- 
ments over the SDO/HMl sequence of magnetograms. The win- 
dow size of 5" is appropriate because it is large enough to con- 
tain the coherent structures to correlate, but not too large to keep 
as much as possible the spatial resolution and also to limit the 
computation time to feasible values. The time series used for this 
analysis starts at 02:00 UT on 17 September 2010 and is com- 
posed by 137 images with a cadence of 45 seconds. Images are 
grouped in 25-min series for the LCT analysis. Prior to applying 
LCT the sequence of images is aligned to eliminate possible jit- 
ter and rotation of the observed target within the field of view. 
Since the differential rotation is of minor influence here due to 
the small field of view, we do not remove this effect. Taking the 
first image as reference, the second image is aligned to it by cor- 
relating them. Then, the third image is aligned with the second 
one with the same method. The process is repeated until all the 
137 images are aligned with each other. We find that the maxi- 
mum shift is up to 21 pixels in the East- West direction and less 
than 1 pixel in the North-South direction. The absolute values 
of the transverse velocities should be taken with caution because 
the LCT technique in general produces so me systematic errors 
in the determination of proper motions (INovember & Simon! 
1988). The LCT method may underestimate 20%- 30% of the 
velocities in extreme cases (lMolownv-Horaslll994b. But it usu- 
ally generates reliable flow patt erns ([November & Simon, 1988; 
Vargas Dommguez et al.Ll2008l) . 

The results confirm the emergence pattern seen in the mag- 
netic field movie. There is indeed a diverging flow with the po- 
larities pi and p2 moving mostly westward and their negative 
counterparts, nl and n2, moving mostly eastward (this diverging 
flow pattern is, in average inclined by about 20° on the East- 
West direction), while the bipoles nl-pl and n2-p2 are mostly 
East- West oriented at that time. 



2.3. Vector Magnetic Field 

We further analyze the vector magnetic fields observed by 
SDO/HMl, which records six flltergrams at each of six wave- 
lengths every 135 s. The flltergrams at each wavelength corre- 
spond six polarization states, i.e., I ± S , where S = Q,U, and 
V, where the notations /, Q, U, and V represent the four Stokes 
parameters. In order to increase the signal to noise ratio and re- 
move the p-mode signals, all the flltergrams are averaged over 
12 minutes. More precisely, the average uses a cosine-apodized 
boxcar with an FWHM of 720 s. The tapered temporal window is 
actually 1215 s. The Stokes parameters, /, Q, U, and V, are com- 
puted from the flltergrams at each of the six wavelengths. The 
vector magnetic field and other thermodynamic al parameters are 
fitted by the inv ersion code of Ver y Fast Inversion of the Stokes 
Vector (VFISV: lBo"rrero et al.ll201 lb . The 180° ambiguity for the 
transverse components of the vector magnetic field has been re- 
solved by the improved version of the minimum energ y method 
dMetcalfl [19941 iMetcalf et ani2^0^lLek7et^l2009h . We ana- 
lyze all the 1 1 vector magnetic fields between 03:00 to 05:00 UT. 
The selected field of view for the analysis is x e [-300", 190"] 
and y e [-560", -315"]. 

To study the electric current well after the recurrent jets 
and to compare with other observations, we analyze two 
more vector magnetic fields, one observed by SDO/HMl at 
09:36 UT and the other one by the Telescope Heliographique 
pour 1' Etude du Magnet isme et des Instabilites Solaires/Multi- 
Raies (THEMIS/MTR: lLopez Ariste et al. 2000; Bo mmier et all 
12007b . THEMIS/MTR is a spectro-polarimeter with higher po- 
larimetry accuracy but lower cadence than that of SDO/HMl. 
THEMIS/MTR is free from systematic instrumental polarization 
(lLopez Ariste et ail l2000h, while the polarim etric accuracy of 
SDO/HMl is about 0.4% dSchouet alj.l2012al) . THEMIS/MTR 
scanned the solar surface from east to west to cover the field 
of view 161" x 104" in the time range of 09:34-10:38 UT. 
The scan step size along east-west direction is about 0.8", the 
pixel sampling along north-south direction is about 0.2". We 
align the magnetic field observed by THEMIS/MTR with that 
by SDO/HMl using a common feature identification method. 
First, a line-of-sight SDO/HMl magnetic field at the middle 
time (10:06 UT) of the THEMIS/MTR observation is inter- 
polated to the THEMIS/MTR spatial resolution. Then, we se- 
lect some common features on both line-of-sight magnetic 
fields and record their coordinates. Finally, the offsets of the 
THEMIS/MTR data referred to the SDO/HMl data are computed 
with the coordinates recorded in the previous step. 

Since the centers of the field of views for the regions of in- 
terest are not close to the center of the solar disk, the projection 
effect must b e removed. We remove t his effect with the method 
proposed by iGarv & Hagvardl (ll990h . which converts the line- 
of-sight (B() and transverse components (Bg and B n ) to the heli- 
ographic components (B x , B y , and B z ) and we project those fields 
from the image plane to the plane tangent to the solar surface at 
the center of the field of view. This coordinate transformation is 
applied in the field of views given above, while a smaller region 
is cut and displayed in Figures|4^ and[4p as an example. 

2.4. Electric Current Evolution 

Using the vector magnetic field, we compute the vertical electric 
current density with the Ampere's law, 

1 IdB, 8B X \ 
yu \dx dy ) 
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Fig. 4. (a) SDO/HMI vertical com- 
ponent of the magnetic field (B z ). 
The projection effect is removed, (b) 
SDO/HMI vector magnetic field with 
the 180° ambiguity being removed. 
(c)-(g) Distribution of the vertical elec- 
tric current density j- averaged over 12 
minutes at different times and observed 
by SDO/HMI. The middle time of the 
observations is indicated at the top of 
each panel. Gray/black lines mark the 
contour levels of the positive/negative 
polarities. The contour levels of B- are 
-1200,-900,-600,300, and 600 G. 
The color bar of j- is shown on the right 
side, (h) Distribution of j- observed by 
THEMIS/MTR centered at 09:42 UT. 
The scan duration is about 16 minutes 
in this field of view. The contour lev- 
els of B : are -882,-661,-440,110, 
and 220 G, which are at the same 
percentage of their minimum (for the 
negative contours) and maximum (for 
the positive contours) values as that in 
Figure |4g. 



where jjQ = 47r x 10~ 3 G m A -1 . The vertical electric cur- 
rent densities are computed for all the analyzed vector magnetic 
fields. Some selected j z distributions before, close to the three 
peaks of, and after the recurrent jets are plotted in Figures |4j;- 
|4ji. The j z distributions close to the jet peaks have very simi- 
lar patterns as shown in Figure |4j;-@f. However, the magnetic 
field at 09:36 UT well after the jets has much lower electric 
current density (Figure |4g) than those in the jets duration do 
(Figure |4j;-f). The j z distribution computed by the vector mag- 
netic field observed by THEMIS/MTR as shown in Figure |4ji 
also shows very low electric currents. The difference between 
Figure |4g and |4ji is caused by the different spatial and tempo- 
ral resolutions and spectro-polarimetric accuracies of SDO/HMI 
and THEMIS/MTR. Besides the above reasons, they are instru- 
ments of two different types, where THEMIS/MTR is a spectro- 
graph and SDO/HMI is a filter/graph. THEMIS/MTR records the 
Stokes line profiles along a slit at one exposure but needs to scan 
perpendicularly to the slit to cover a required field of view. While 
SDO/HMI records an image at a wavelength point at one expo- 
sure but needs to scan along the wavelength to get the Stokes 
line profiles. 

Next, we integrate \j z \ in the field of view as shown in 
Figures |4^-0Ji for those regions where \j z \ is above the noise 
level 0.02 A m~ 2 estim ated by the formula (1 /fj.o)(5Bj/ Ax) 
dGarv & Demouhnl. 11995b . where 5B T ~ 100 G is the assumed 



transverse field error and Ajc = 0.5" is the grid size. The field of 
view in Figures @k-|4ji are selected to be similar to the one used 
to compute the line-of-sight magnetic flux (Figures QJl and|2j)). 
But they are not identical, since the projection effect has been 
corrected to compute the vertical electric currents. Figure^ dis- 
plays the evolution of the integrated vertical electric current, 
where an obvious quasi-periodicity can be found. This recurrent 
evolution is present on the background of a stationary current 
pattern, whose average electric current is (2.51 +0.06) x 10 12 A. 
The current pattern is stationary since the vertical current den- 
sity distribution does not change too much during the jets as 
shown in Figures |4j;-[4f . The fluctuation of L around the back- 
ground (defined as |max(/ z ) - min(/,)|// z , where /, denotes the 
average of I z ) is limited within 8%, which is above the errors 
of I z (around 2%). The integrated electric currents, I z , for the 
magnetic fields observed by SDO/HMI at 09:36 UT (Figuregg) 
and by THEMIS/MTR at 09:42 UT (Figure gji), are around 
(1.88 + 0.05) x 10 12 A and (0.43 + 0.04) x 10 12 A, respectively. 
Therefore, the currents significantly decrease later on. 

The errors for the integrated L as shown in Figure^ are es- 
timated by a Monte-Carlo method. The basic assumption is that 
the errors for I z are normally distributed in the whole field of 
view of Figure |4] and three times the standard deviation of the 
errors is 0.02 A m -2 . We add those errors to I z and integrate 
their absolute values with the criteria described in the previous 
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paragraph. After repeating 50 times for each observation time, 
the errors are estimated as the standard deviation of the 50 inte- 
grated 

Figure [5] also shows the integrated vertical electric currents 
overlaid on the 171 A flux. The integrated vertical current is de- 
rived from the de-projected heliographic vector maps, while the 
171 A flux is computed in the image plane (plane of the sky). 
The de-projection allows to derive the vertical current density, 
which has a more physical meaning than the line-of-sight one 
before de-projection. Since the EUV jets are coronal events and 
the electric current is on the photosphere, the field of view for 
computing the vertical electric current does not need to be the 
same as that for the 171 A flux. We find that the peaks on the 
I z curve have a good coincidence with those peaks on the 171 A 
flux. They both repeat three times. However, the magnitudes of 
the peaks on the /, curve have large differences, which can be ex- 
plained by the lower temporal cadence of SDO/HMl vector mag- 
netic field observations (12 minutes). Therefore, sample points 
on the I z curve are randomly distributed compared to the 171 A 
flux peak time (Figure^). 

In order to mimic a higher temporal resolution of the I z evo- 
lution, we combine the observations at the three peaks with one 
assumption. We assume that the vertical electric current accu- 
mulation and decrease processes are similar for all the three re- 
current jets, which are homologous from EUV observations, i.e., 
we assume that the peak magnitude of the I z curve and the evolu- 
tion timing relative to the 171 A flux are almost the same. Then, 
we shift the I, curve around t\ = 03:17 and ?3 = 04:22 UT to the 
time ?2 = 03:57 UT as shown in Figure^. From the combined 
I z curve in Figure \5fa, we find that the vertical electric current 
increases first and then decreases in a jet process. 

The combined I z curve in Figure [5j> also shows two addi- 
tional features. First, the I z curve peak is earlier than the 171 A 
flux peak. Second, the vertical electric current seems to increase 
faster than it decreases. However, we remind the limited tempo- 
ral cadence of SDO/HMI, so both results will need to be checked 
with other observations. 



2.5. Magnetic Configuration before the Recurrent Jets 

The magnetic configuration is one of the necessary informa- 
tion for discriminating different magnetic reconnection models. 
For this reason, we study the three-dimensional magnetic field 
with the nonlinear force-free field (NLFFF) extrapolation of the 
coronal field from the vector magnetic field as the boundary 
condition. The optimization method is selected for the N LFFF 
extrapolation (IWheatland et all l2000t IWiegelmaii nl 12004 . The 
bottom boundary is provided by the SDO/HMl vector mag- 
netic field at 03:00 UT, so just before the studied jets. The pro- 
jection effect has been corrected as described in Section 12.21 
Since the photospheric vector magnetic field does not satisfy the 
force-free condition, we apply the preprocessing method to re- 
move the net magnetic fo rce and torque on the bottom boundary 
dWiegelmann et all 120061) . A requirement for the preprocessing 
is that, the magnetic field at the bottom boundary needs to be 
isolated and flux balanced. Therefore, we select a field of view 
as shown by the dotted box in Figure flj to enclose the major 
magnetic flux of AR 1 1 106. Since removing the projection effect 
changes the geometry of the field of view, we cut off the edges 
where the data are incomplete to get a rectangle box, which is 
resolved by 932 x 486 grid points with Ax — Ay « 0.5". In or- 
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Fig. 5. (a) Evolutions of the 171 A flux (normalized as in 
Figure |2]i and the unsigned vertical electric current I z , which is 
integrated in the field of view as shown in Figure |4] Blue stars 
mark the normalized 171 A flux and red diamonds mark the un- 
signed vertical electric current /,. (b) The electric current evo- 
lution curves at the first and third peaks have been shifted to 
the second peak and the total profile is shown with a continuous 
line. The shifted time differences are ?2 - 1\ and tj, - ?2 for the two 
curves (denoted by dashed and dash-dotted lines), respectively. 
Note that the vertical axes for both the 171 A flux and I, do not 
start from zero. 



der to measure the flux balance quantitatively, we define the flux 
balance parameter as 



Zi \(B z )i\dS ' 



(2) 



where the pixel index i runs all over the field of view, and dS is 
the area of one pixel. The flux balance parameter e/ is about 0.03 
for the magnetic field in this field of view, which is a tolerably 
small imbalance. 

The boundary field of view for the NLFFF extrapolation is 
shown in Figure^. There are two reasons why we use this small 
field of view. First, the region of interest where the jets originate 
is small. It is acceptable if we put it in the center. Secondly, the 
spatial sampling of SDO/HMl is high (0.5"). Thus, in the se- 
lected field of view (about 100" x 100") and the height range, the 
computation box includes 200 x 200 x 200 grid poi nts, which are 
reason able sizes for the computation. However, DeRosa et al. 
(l2009h pointed out that it is critical for the vector magnetic field 
to cover as large areas as possible for successful NLFFF mod- 
eling. In order to fulfill this requirement in the limited field of 
view, we use the following way to prepare the boundary and ini- 
tial conditions for the higher-resolution NLFFF extrapolation. 
First, we compute a lower-resolution NLFFF with the vector 
magnetic field in the field of view used above for the prepro- 
cessing as shown by the dotted box in FigureQJ. Due to the lim- 
itation of computation resources, the lower-resolution NLFFF 
is computed in a box of 466 x 243 x 201 grid points with 
Ax — Ay = Az ~ 1.0". Then, we cut out a sub-volume in the field 
of view as shown in Figure [6^ with a height of 100" and interpo- 
late the NLFFF to the spatial resolution of Ax = Ay — Az ~ 0.5". 
The interpolated NLFFF is used as the initial, lateral, and top 
boundary conditions for the higher-resolution NLFFF extrapo- 
lation. The bottom boundary is provided by the aforementioned 
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X (arcsecs) 



Fig. 6. (a) Magnetic field lines computed from the nonlinear force-free field model. Different colors indicate different field line 
systems. The background is the vertical magnetic field B, shown with grey levels. White/black color represents positive/negative 
magnetic polarity. The box marks the field of view shown in panels (b)-(f). (b, c) Orange ribbons represent the QSL sections on the 
photosphere. Only the QSLs with the squashing degree Q > 10 14 are plotted. Grey/black contours are the same as that in Figure Et. 
(d, e) Magnetic field lines overlaid on the photospheric distribution of the vertical electric current density j z as shown in Figure 2}:. 
The current j z is computed with the non-preprocessed vector magnetic field, (f) Similar to panels (d) and (e) but only some sample 
field lines with magnetic dips are shown. 



preprocessed vector magnetic field in the field of view as shown 
in Figure [6^, which is cut from the larger flux-balanced area. 
We use a buffer zone with 25 grid points, in the NLFFF opti- 
mization method, to decrease the effect of these boundaries. The 
weighting function decreases from 1 to with a cosine profi le in 
the buffer zone. Fin ally, we refer to IWheatland etall (l2000t> and 
Wiegelman iil (l200l for more detailed descriptions of the extrap- 
olation method. 

We adopt two metrics to test if the NLFFF has reached the 
force-free and divergence-free state. The force -free metric is de- 
fined as the current-weighted average of sin (IWheatland et al.L 
2000): 



< CW sin 6 >-- 



where 



sin 0j 



IJ/ x B ; 

JiBi 



(3) 



(4) 



represents the angle between the current density J and the mag- 
netic field B, and B = |B|, J = |J|. The summation is done in 
the sub-domain of 77 x 75 x 35 grid points, whose projection is 
shown as the solid box in Figure [6^. Only the grid points where 
/ > 0.02 A m -2 (the noise level in Section |2~4] | are considered 



in computing sin 6. For a perfectly force-free magnetic field, the 
force-free metric should be zero, i.e., the current density J and 
the magnetic field B are parallel to each other. The divergence- 
free metric is defined as the unsigned average of the fractional 
flux unbalance, < >, where 



\fi\ 



l(V ■ B),-| 

6Bi/Ax 



(5) 



The divergence-free metric is computed in the same domain as 
that for the force-free metric, while the difference is that all the 
grid points of the sub-domain are considered. This metric should 
be zero in a truly divergence-free field. 

The force-free and divergence-free metrics for the NLFFF 
model are 0.19 and 3.0 x 10~ 3 , respectively. The correspond- 
ing angle for the force- free metric, which is defined as arc sin < 
CWsin6> >, is 11.0°. iMetcalf et al.1 J2008) compared various 
NLFFF algorithms using simulated chromospheric and photo- 
spheric vector fields. The force-free metric 0.19 that we de- 
rive here is slightly better than that of 0.26 using the prepro- 
cessed and smoothed photosphe ric boundary with the optimiza- 
tion code o fWiegelmann(2004). However, it is larger, therefore 
worse, than the force-free metric 0.11 using the chromospheric 
boundary. The divergence-free metric 3.0 x 10~ 3 is larger (and 
worse) than both the metrics using the chromospheric and pre- 
processed and smoothed photospheric boundaries with the op- 
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timization method as listed in Table 5 of Met calf et alJ (l2008h . 
since they adopted simulated vector magnetic fields as boundary 
conditions, which are smoother than the observational bound- 
aries that we use in this paper. 

Some selected magnetic field lines in the NLFFF model 
are plotted in Figure [6] The overall magnetic configuration in 
Figure [6^ shows different connectivities with small closed field 
lines linking the small diverging bipole (represented by the red 
lines connecting positive polarities pi and p2 to negative ones 
nl and n2) and long field lines (represented by the blue lines an- 
chored in polarities NO and Nl). We find no magnetic null point 
in between these two set of field lines, so no separatrix. Rather, 
a continuous but drastic change of magnetic connectivity, so a 
QSL, is separating both types of c onnectivities. We compute the 
squashing degree Q as defined in lTitov et al.l(l2002l) to locate the 
position of the QSLs. Figures |6p and [(J show the QSL sections 
on the bottom boundary. The closed (red) and long (blue) field 
lines are clearly separated by QSLs between them. 

We also overlay the magnetic field lines on the photospheric 
distribution of the vertical electric current density j, as shown 
in Figures [6}l and|6^. The footpoints of those field lines close to 
the QSL are rooted on the regions where j z is large. This result 
is coherent with previous findings where concentrated electric 
currents have b een found at the border of the QSLs in flare con- 
figurations (see Demou lin et all IT9971 and references therein). 
For a thin QSL, a moderate difference of magnetic stress across 
the QSL creates strong electric current. Then, a current layer is 
easily build up at a QSL. For the first time, the high temporal ca- 
dence and magnetic sensitivity of SDO/HMl allows to diagnose 
such current layer evolution for such small events as the studied 
jets. 

Due to the existence of the parasitic positive magnetic po- 
larities, some blue field lines are bent down from higher al- 
titudes to lower altitudes. Magnetic dips and bald patches are 
present where the magnetic field lines are bent down as shown in 
Figure |6f. We refer to the following theoretical or observational 
papers for detailed de finition and analys is of magnetic dips an d 
bald patches, such as iTitoy et al.1 (1 1993b. iBungev et alJ d 1996b . 
Man drini et al.l (12002b . and lDemourinl(l2005b . However, it is not 
clear if the bald patches played any role in the magnetic recon- 
nection of the studied jets, which asks for further studies. 

3. Discussion and Conclusion 

In summary, we find that three EUV jets observed by SDO/AIA 
from 03:00 to 05:00 UT on 2010 September 17 recurred in the 
same region close to the border of AR 11106. The time differ- 
ences between the three successive 171 A peaks are 40 minutes 
and 25 minutes. On the border of the following negative polarity, 
parasitic positive polarities (polarities pi and p2 in Figure [3} are 
present close to the footpoints of the EUV jets. They are parts 
of a magnetic bipole which is growing in size as shown by the 
diverging flow pattern. The line-of-sight magnetic flux evolution 
does not have a clear relationship with the 171 A flux evolution. 
However, the high time cadence of SDO/HMl (12 minutes for 
the vector magnetic field) allows us to follow the photospheric 
current evolution. We integrate the absolute vertical electric cur- 
rent density found in the region of the jets. The peaks of this total 
current have a good temporal coincidence with the peaks of the 
171 A flux. 

The aforementioned feat ures can only be par tly explained 
by the emerging flu x model dShibata et all 1992b) and the con- 
verging flux model dPriest et al.lll994|) since the newly emerged 



magnetic flux is consistent with the former model and the 
bald patch configuration is consistent with the latter one. But 
both models are two dimensional with magnetic reconnec- 
tion in separatrices (as implied by their dimensionality), while 
the magnetic connectivity is not necessarily discontinuous in 
the three dimensional space. The above models can be gen- 
eralized to thre e dimensional configuration having a magnetic 
null point (e.g..lMoreno-Insertis et all l2008t iTorok etaQ [2009 ; 
iPariat et all [2009. 2010). However, in the present study, no null 
point is found in the corona with the NLFFF extrapolation. 
Rather continuous changes of connectivity are present with dras- 
tic changes at some locations. Indeed, the computed coronal con- 
figuration has a QSL separating a small-scale evolving double 
bipole (polarities pl-nl and p2-n2 in Figure [3} from the large 
scale AR magnetic field, a configuration comp arable to those 
found previously, e.g. in an X-ray bright point (Man drini et all 
[l996h . 

The temporal evolution of the photospheric magnetic field 
(see the movie attached to Figure [3) as well as the velocities 
deduced by LCT both indicate divergent flows in the small par- 
asitic bipoles. These flows would force the closed field lines of 
the small parasitic bipole to grow in size, and to interact with the 
overlaying large scale field lines (Figure |6). Such kind of evolv- 
ing magnetic configuration is known to build up a narrow current 
layer in the QSL as conjectured analytically (Dem oulin et all 
1996) and found in n umerical simulations even for relatively 
small displacements (lAulanier et all 120051: lEffenberger et all 
2011). The main reason is that the magnetic stress of very dis- 
tant regions, generated by photospheric flows, are brought close 
to one another, typically over the QSL thickness. This current 
layer becomes thinner and stronger with time. When the currents 
reach the dissipative scale, there is a breakdown of ideal MHD 
in the QSLs, magnetic re connection occurs a nd part of the mag- 
netic energy is released (lAulanier et alll2006l) . In the above stud- 
ied configuration, this process would create some newly formed 
lower-lying field lines with less shear and some newly formed 
higher field lines, along which hot collimated plasma are ejected 
into the higher corona. 

The quasi-periodicity and homology of the coronal jets re- 
quires two necessary conditions. First, the magnetic energy is in- 
jected into the corona uninterruptedly in almost the same region, 
and secondly, it releases intermittently. The first condition is sat- 
isfied by the observed diverging photospheric motions which in- 
ject free magnetic energy in the coronal field without changing 
too much the photospheric magnetic field distribution. The sec- 
ond condition requires the current layer to become thin enough 
recurrently in order to reconnect and release rapidly part o f the 
plasma trapped in the low lying loops (of the diverging bipole). 
The observed periodic build up and decrease of the photospheric 
currents, in phase with the EUV jets (Figure [5), is an evidence 
that this second condition is met. Indeed, only a small fraction of 
the current is dissipated, an indication that the reconnection stops 
rapidly when the current layer weakens. Then, the next build up 
of current starts from a non-potential configuration. This allows 
a relatively fast build up of a new thin current layer and may 
explain the small time interval between the jets (40 and 25 min- 
utes). The photospheric flows displace the magnetic polarities 
by only about 300-480 km during the successive jets, where we 
use v = 0.2 km s _1 estimated by the LCT velocity data as the 
average velocity in the encircled region of Figure [3] 

What is still surprising is that the build-up and decay times 
of the currents are comparable (Figure |5). Usually, the build- 
up time is thought to be much longer than the decay time for 
many solar activities, such as flares and coronal mass ejections 
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(CMEs). A reason for the comparable time scale could be the 
limited time cadence of SDO/HMl including the constrain of a 
large enough signal to noise ratio, so this kind of study needs a 
higher cadence and better polarimetric accuracy to get a more 
precise temporal evolution of the currents. We mimic a higher 
temporal cadence by supposing that the evolution processes are 
the same, then we combine the current measurements of the 
three jets using the EUV flux maximum as a time reference. 
With these limits, we find a current maximum near the begin- 
ning of the EUV flux rise, about ten minutes before the EUV flux 
maximum. This indicates that the free magnetic energy starts to 
decrease as soon as there is an evidence of reconnection. 

In conclusion, we find evidence indicating that the studied 
recurrent coronal jets are caused by magnetic reconnection in the 
QSL present between a diverging bipole and the main AR mag- 
netic field. Magnetic energy is injected by the continuous diverg- 
ing motions observed at the photospheric level. We assume that 
this process has build a thin current layer whose photospheric 
cross sections is deduced from the vector magnetograms. The 
current evolution is found in phase with the EUV flux with a cur- 
rent maximum present before the jets followed by a relaxation 
to lower value. However, this evolution of the current intensity 
is limited to around 8% of the average background. Note that 
only the vertical component of electric currents are derived on 
the photosphere. If the NLFFF assumption is valid, the coronal 
currents are directly related to the photospheric currents , since 
a = fioj z /B z and a is constant along a field line. Therefore, we 
conclude that the magnetic system always stays close to the re- 
sistive instability of the current layer. The jets does not change 
the photospheric vertical current distribution too much and only 
a moderate photospheric evolution is needed to restart a new jet. 

All together, we interpret these observations as the conse- 
quence of the continuous photospheric driving of the coronal 
field. A thin current layer is built in the QSL, where the magnetic 
field reconnects when the electric currents are strong enough and 
stops reconnecting soon later on. The magnetic reconnection dis- 
sipates only a small part of the electric currents. The above phe- 
nomena repeat in time and drive the recurrent jets as the confined 
plasma in the closed loops is accelerated, after reconnection, into 
the large scale field. 
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Fig. 7. Online only figure. A movie showing the evolution of this 
figure from 03:00 to 05:00 UT is available in the online journal, 
(a) SDO/AIA 171 A image. The box marks the region where the 
171 A fluxes are computed, (b) SDO/AIA 171 A flux normal- 
ized to the background level, which is defined as the flux in the 
quiet Sun region (see Section l2?T1 for more details). The vertical 
dashed line marks the time of the 171 A image in panel (a). 
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